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^ : 

^""1 . The observable Universe is described by a collection of equal mass galaxies linked 

C " ^ 1 into a common unit by their mutual gravitational interaction. The partition func- 

' tion of this system is cast in terms of Ising model spin variables and maps exactly 

C " 3 ' onto a three-dimensional stochastic scalar classical field theory. The full machin- 

^ ery of the renormalization group and critical phenomena is brought to bear on this 

^\ ■ field theory allowing one to calculate the galaxy-to-galaxy correlation function, 

■""^ 1 whose critical exponent is predicted to be between 1.530 to 1.862, compared to the 

phenomenological value of 1.6 to 1.8 

u'. 

bJQi It is an empirical fact that the Universe today is dominated by matter 

^ I which is in gravitational interaction and is organized into spatial configura- 

" '~j ■ tions consijating of galaxies, groups of galaxies, superclusters, and even larger 

rS I structures Elia. It is also true that one of the outstanding problems in mod- 

^ . ern cosmology is to provide a quantitative description of this structure from 

■ ■ ■ ' first principles u. One of the standard ways in which structure, that is, the 

departure from randomness and homogeneity, is described, is by means of cor- 
relation functions. These objects are a basic tool for unravelling structure in 
physical systems ranging from gases and liquids to solids. Experimental input 
is typically deduced from analysis of the scattering of neutrons or x-rays. By 
contrast, in cosmology, this input is provided by direct observation. From ob- 
servation, the galaxy-to-galaxy correlation function (Cgaiif)) is known to scale 
as ^(r) ~ r~'' , wherejy ~ 0(1.6 to 1.8), instead of 7 = 1, which is what one 
would naively expect □ for a homogeneous distribution of matter in a three 



dimensional space. 



Based on a talk presented by D. Hochberg at the International Workshop on Elementary 
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Given these observational facts, we shall demonstrate how to calculate the 
exponent 7 from first principles U. We take our cueJpy analogy with the well 
established theory of liquids and interacting gases where the relative ar- 
rangements of the atoms and molecules (galaxies, in our case!) are accurately 
described in terms of the principles and methods of classical statistical me- 
chanics. We will apply the same techniques to a system made up by many 
galaxies in gravitational interaction, subject to fluctuations emerging from the 
intrinsic properties of the gravitational interactiond. Wc will obtain the canon- 
ical partition function for this system and compute the two-point function and 
its corrections due to these fluctuations. 

Begin by considering the continuous mass density p{r) describing the 
spatial distribution of galaxies. The departure from perfect homogeneity is 
Sp{r) — p{r) — p, where p denotes the average density. The galaxy-galaxy 
correlation function is definedlj as ^Gai(rj — =< Sp{ri)dp{rj) > /p^, the 
angular brackets denoting a suitable averaging. This is the joint probability 
distribution function for finding a galaxy at position given that there is one 
at Tj. In the static, nonrelativistic and weak gravitational field limit, (we dis- 
cuss the effects of expansion later) the gravitational interaction energy for this 
system is 

H^nt = J y"'^'^i-,d3rjp(r,)^^:--^p(rj), (1) 

where the integral excludes the point = Vj, and the factor of 1/2 appears 
because each galaxy is counted twice in the integration. It is reasonable to 
consider the gas of N galaxies as made up of discrete spatially localized points 
of equal mass nio at position r^, so that p{r) = J2f=i''^jH^ ~ ^j)- The 
interaction energy then becomes 

where rrij is the mass of the galaxy at site r^, and i ^ j- In this "granular" 
picture, the contrast function 5p{vi)/ p is equal to 1 if there is a galaxy at r,; and 
is equal to —1 if there is a void. Because we have assumed all N galaxies have 
equal mass, there is a natural correspondence between the irii and a two- valued 
(±) "spin" variable 

mi = ^{si + l). (3) 

A galaxy at site r; corresponds to "spin-up" , while a void or deficit corresponds 
to "spin-down". Because of (3), we can express the Boltzmann factor and 
associated partition function of this gravitational Ising system in terms of the 
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spin variables as: 



where A^^^ — , and hj — \ J^-^\ -i ^'^'^ have dropped an irrelevent 

constant term. /3 plays the role of an inverse temperature. The sum in (4) is 
over all possible spin configurations si — ±1, S2 = ±1, ■ ■ ■ , sn = ±1- 

To continue, we make use of the Hubbard-Stratonovich (or gaussian) trans- 
formationQ 

{detA)-^/^ exp 



(5) 



\ m,n m J 

which maps exactly the discrete system described by the spin variables Si onto 
a system described by the continuous real variable — oo < Xm < oo. Applying 
(5) to the partition fimction in (4) yields 

ZfsZim ^ Cj m exp[-^ J (0(r) - hir))A{r, r')(0(r') - h{r'))drdr' 
rfr logcosh(/30(r))] (6) 



performing the exact configurational sum over the {si} and after taking the 
continuum limit which rigorously converts the discrete system into a continu- 
ous three dimensional 1-component scalar field theory. The linear term in (4) 
results in the highly nonlinear log cosh appearing in (6). C is an inessential 
factor and h{r) = ^ J dr' A~^{r,r'). The continuous scalar field (j) is the order 
parameter for this system; its n-point Green functions correspond to the sta- 
tistically averaged product of n times the contrast evaluated at n-sites. The 
operator A in the continuum limit is given by 

A(r,r') = <5(3)(r-r')— (7) 
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since = -inS^^Hx). 

As is well knownH, the connected two-point function for the spin system 
and the field theory are the same. Moreover, the unavoidable fluctuations in 
the field result in an anomalous dimension rj shifting the canonical dimen- 
sion of the field such that when |r — r'| oo, the two-point function for the 
hamiltonian in (6) scales as 

lim <s,s, lim ecaidr-r'l) ^ |r-r'r*''"'+"^ (8) 

|r— r'l — ^oo |r — r'l — ^oo 

Here d is the dimension of space {d = 3) and r/ is the critical exponent for 
the pair correlation fimction whose value (0.0198 — 0.064) (see Table 1) differs 
from zero due to the fluctuations in 0. 

Thus, our calculation shows that for large separations, the galaxy-galaxy 
correlation function must scale as 

^Gai{h-r,\)^r-^''-^+^\ (9) 

withl 1.0198 < 1 -I- 77 < 1.064 for d = 3, and without taking the background 
expansion into account. 

Thus far, our calculation has been performed for a static ensemble of galax- 
ies, but the Universe is expanding and the effects of this expansion will modify 
the values of the critical exponents. From the theory of dynamical critical 
phenomena, it is phenomenologically known from computer simulations that 
the way time enters into the correlation function is by a modiflcation of its 
argument, viz., 

ar)^^{r;t)^Fir/lit)), (10) 

where l{t) ^ t'' is a. function of time and here, r denotes the comoving coordi- 
nate, r comoving 1 which is related to the physical coordinate via 

'^physical ^{^^ ^comoving : (-^-^) 

with a{t) the scale factor for the background spacetime. From computer simu- 
lations in condensed matter physics, it is established that ^ = 1/3 for systems 
with a conserved order parameter, while ( = 1/2 when the order parameter is 
not conserved. 

For a matter dominated expansion, the case we are interested in, a{t) = 
t^/^ and thus for a given physical separation 

t ~ {^T'comoving) ^ • (-^^) 

This means that the function l{t) scales as 

^(^) — — ij" coraoving^ • (-^^) 
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Table 1: Calculated values for the galaxy-galaxy correlation function critical exponent 



Method of Calculation 


IStatic 


^Expandinq 


1 Expanding 


Scries estimates 


1.056 ±0.008 


1.584 ±0.012 


1.848 ±0.014 


0{e) 





1.5 


1.75 


1.0198 


1.530 


1.785 




1.037 


1.555 


1.815 




1.029 


1.543 


1.801 



Putting together these facts yields the scahng behavior of (,{r;t): 



^(r, — I T(2onioving / 



-3C/2 



comoving/ ' comoving 
comoving 



-(d-2+ri) 



(r . )-{d-2+ri){l+3<:/2) ^i.N 

y comoving J ■ 



In other words, the critical exponent resulting from the expansion of the Uni- 
verse and dynamical critical phenomena effects is 

7= (d-2 + r;) x (1 + 3C/2). (15) 

The predicted (calculated) value for 7 is between 1.530 and 1.596 for ( — 1/3 
(for a conserved order parameter) and between 1.75 to 1.862 for C = 1/2 (non- 
conserved order parameter). Table 1 displays these exponents for both the 
static and expanding cases (C, NC, denote conserved and nonconserved order 
parameter, respectively). The static values are taken from Ref[9]. The last 
four entries refer to expansions in powers of e = 4 — d. These values are to be 
compared with the values inferred from current galaxy iatalogs, which range 
from 1.5 for the APM survey, to 1.8 for the Lick surveyHu. 

We have thus accomplished what we set out to do, namely, provide a first- 
principles calculation of the galaxy-galaxy correlation function using standard 
methods of statistical mechanics and the theory of critical phenomena. Nat- 
urally, one might (and should) ask why this calculation is successful, in other 
words, what does a system of galaxies have to do with a spin model? The 
answer is provided by the universality hypothesis, according to which rather 
distinct physical systems will exhibit identical critical behavior provided the 
systems in question possess the same space dimensionality and have order 
parameters with the same dimensionality (or number of independent compo- 
nents). We see that this is the case here: d — 3 for both systems, while the 
scalar field and the magnetization are one-component order parameters for 
the galaxy system and the 3d Ising model, respectively. In other words, both 
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systems belong to the same universality class. Thus, we must expect them to 
have the same values for their critical exponents. 
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